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Abstract. Lie group theory states that knowledge of a m-parameters 
solvable group of symmetries of a system of ordinary differential equa- 
tions allows to reduce by m the number of equations. We apply this 
principle by finding some affine derivations that induces expanded Lie 
point symmetries of considered system. By rewriting original problem in 
an invariant coordinates set for these symmetries, we reduce the number 
of involved parameters. We present an algorithm based on this stand- 
point whose arithmetic complexity is quasi-polynomial in input's size. 



1 Introduction 

Before analysing a biological model described by an algebraic system, it is useful 
to reduce the number of relevant parameters that determine the dynamics. 

Example 1. In order to give an example of such a reduction, let us consider the 
following Verhulst's logistic growth model with linear predation (see § 1.1 in [1]): 

x = (a — bx)x — cx, a = b = c = 0, t = 1, (1) 

for which all forthcoming computations could be easily performed by hand. As- 
suming that a ^ c, one can represent the flow (t, x) of (1) using parameterization: 

t = t/{a-c), x = (a-c)x/b, (2) 

where (t,x) is the flow of the following differential equation: 

x=(l-x)x. (3) 

In this formulation of (1), parameters a and c were lumped together into o — c 
and its state variables x and t were nondimensionalise. 

Usually, presentation of this kind of simplification relies on rules of thumbs (for 
example, the knowledge of units in which is expressed the problem when dimen- 
sional analysis is used) and thus, there is — up to our knowledge — no complexity 
results on these kind of reduction methods (see [2] and references therein) . 

However, these reductions are generally based on the existence of Lie point 
symmetries of the considered problem (for reduction based on dimensional anal- 
ysis, see the Theorem 3.22 in [3]). 



Example 1 (continued). The following continuous groups of transformations: 



a — A, i a — ► V a , 

S M : ' < -» (4) 

\ K UK., 

c — A, c — > Vc, 

leave invariant system (1) and its solutions. These symmetries are called ex- 
panded because they act on the expanded space of variables that includes the 
system parameters in addition to independent and dependent variables. 

The system (3) is obtained by 'factoring out' the actions of the symmetries (4) 
of the original system (1) and thus, it is invariant under these actions. The 
relations (2) parameterize the solutions of the original system (1) in function of 
the solution (t,x) of invariant system (3) and of the free parameters a, b, c; they 
are defined by the composition (Sr„/ (i „\ o X)(t,x, a, b, c). 

The aim of this note is to show how Lie theory unifies and extends the classical 
methods (exact lumping, dimensional analysis, etc.) used to simplify algebraic 
(differential) systems. We adopt a presentation based on algebraic tools closer 
to the actual computations (mainly Jordan normal form and linear algebra) on 
which are based our reduction process. 

1.1 Related works 

The literature on investigation of the invariants of Lie group action and their 
applications to algebraic systems is far too vast to be reviewed properly here. 
The book [4] shows various applications of invariants in the study of dynamical 
systems under a computer algebra viewpoint. Section 4.1 for example, shows how 
the knowledge of invariants of a given dynamical system could simplify further 
computation on it by reducing the degree of involved polynomial expressions. 
In [5] authors show, given a rational group action, how to compute a complete 
set of its invariants. 

We adopt the same general philosophy — determine some system's symme- 
tries and use their invariants — but our purposes is more to reduce the number 
of variables involved in these expressions than their degrees (that is for us only 
a byproduct). Furthermore, while the symmetries considered in [4,5] are quite 
general — and thus, their computation of invariants are exclusively done using 
Grobner bases computations — we restrict oursclf to the use of affine Lie sym- 
metries and thus, the required operations are restricted to linear algebra over a 
number field and univariate polynomial factorization. 

1.2 Main Steps and Tools of the Reduction Process 

The reduction process introduced in this note is based on classical Lie theory but, 
to the best of our knowledge, it is not described elsewhere in the literature. Let 
us present now the tools used in this process through our introducing example: 

Example 1 (continued). Reduction of Example 1 is done as follow 




Step 1. Determine affine infinitesimal generators that induce expanded Lie symme- 
tries of the considered system. If there is no such derivations, our reduction 
process stops. In our example, they are: 

Si--=k + h> 6 2--= x m- b m> b^afc + bfa + cfc-tfc. (5) 

The section 3.1 describes the required computations and shows that their 
complexity is quasi-polynomial in input size. 

Step 2. Choose a generator that is a symmetry of the others (above infinitesimal 
generators form a solvable Lie algebra and our reduction process have to 
take this property into account). As [Si, S2] = [82, £3] = and [Si, 83] = Si, we 
could choose <5i or 8 2 - Determine a principal element (p := log b) of this last 
generator i.e. an element defining a coordinates set in which the derivation 82 
is rectified (equal to the translation d/dp). 

Step 3. The chosen principal element induces an invar iantization of the considered 
system i.e. the system x = (a — x)x — cx that is invariant under the ac- 
tion : — > px. , i — ► i of the one-parameter group of symmetry induced 
by 82- The solution x of the original system is then parameterized by Sb (x) 
where x is a solution of its invariantization and b is a free parameter. 
Repeat Step 1 (supplementary affine symmetries could appear after Step 3). 

Let us stress that it is generally hard to find a general infinitesimal generator 
of a system's symmetry (Step 1) and to give an explicit representation of an 
invariant coordinates set (Step 2) for it. Thus, we restrict ourself to Lie symme- 
tries associated to affine infinitesimal generators for which invariant coordinates 
computation is easy (for general case see [5] and references therein). Hence, we 
do not follows methods developed for general cases because their complexity 
are likely exponential in input's size while we focus our attention to method of 
quasi-polynomial complexity. 

To conclude, remark that this reduction process works also for purely alge- 
braic system (describing fixed point of a dynamical system for example). 

Outline. In the next section, we recall some basic definitions concerning con- 
sidered systems and related derivations. Then, we present the notion of principal 
element and show how it could be used in order to define a rectifying coordinates 
set for general derivations. In the second part of this note, we focus our atten- 
tion on affine Lie point symmetries in order to propose a probabilistic strategy to 
compute them and their associated principal elements. We show how previously 
introduced notions are used in the reduction process by considering invarianti- 
zation of purely algebraic (resp. differential) system and their parameterization. 
Finally, in conclusion we make some remarks and suggest possible further works. 

2 Considered systems and associated derivations 

2.1 Some Algebraic Systems used in Analysis of Biological Model 

Note 2. Notations — Hereafter, we consider an explicit algebraic ordinary dif- 
ferential system S bearing on n state variables X := (xi, . . . , x n ) and depending 



on I parameters := (61, . . . , 9t): 



X = F(t,X,0), 

i = l, = 0. K ' 

Denoting the set {1, . . . , n} by N, the letter X stands for first order derivatives 
of state variables (ij \ j G N) w.r.t. time t and F := (f 3 \ j G N) is a finite subset 
of K(t, X, 0) where K is a subfield (Q for example) of C. In order to determine 
the qualitative properties of the dynamical system (6), it is usual to consider the 
following systems for various subset J of N: 

x 3 = f 3 {t,X,0), VjGJ, 

f l (t,X,0)=O, VieN\J, (7) 
i = l, = 0, 

in which some state variables are considered as parameters. In fact, for i in N, 
the system E^y defines the so-called a^-nullcline of S and the purely algebraic 
system defines its fixed points (see examples of applications in [1]). 

Remark 3. In the sequel, we are going to avoid — as much as possible — any dis- 
tinction between time, state variables and parameters i.e. we work in an ex- 
panded state space (see [6] for another application of this standpoint); hence, let 
us denote the set (t, X, 0) by Z := (z l \ 1 < i < 1 + £ + n) and its cardinal by m. 

2.2 Infinitesimal Generators, associated flows and their rectification 

First let us recall some basic facts about derivations. 

Definition 4. Given a polynomial algebra a derivation ofK[Z] with con- 

stant field K is an additive mapping 5 : K[Z] — > K[Z] that satisfies Leibniz rules: 

V(fi,fi) G K[Z} 2 , S(hh) = h5h+h5h, (8) 

and have K in its kernel. We denote by T)er^K[Z] the set of all such derivations. 
The Lie bracket is defined by the K-bilinear map: 

[,}: Der K IK[Z] x Dcr K K[Z] Der K K[Z], (9) 
{Si, 62) — ► 5i^2 - S 2 Si- 

This map is skew- symmetric and satisfies the following Jacobi identity: 

V (61,62,63) c Der K K[Z], [61, [6 2 , 6 3 }} + [6 2 , [63, Si}} + [63, [6 2 , Si}} = 0. (10) 

The set Der^K[Z} is a IK vector-space spanned by the set of canonical deriva- 
tions {d/dzi, . . . , d/dz m }. It is also a Lie algebra with Lie bracket as product. 

Remark 5. An algebraic system Sj defined by (7) could be seen as a derivation: 

D J '■= m +T, J ejfo§^' ( n ) 
associated to the algebraic relations {/, = 0, \/i G \ J}. 



Derivations considered as infinitesimal generators. The exponentiation 
of a derivation 5 induces several morphisms as shown by following definitions: 

Definition 6. Given a derivation S and r one of its constant (St = 0), one 
can define the exponential map e rS :— X^eN T% ^ % l % ^ from M.[Z] into the alge- 
bra K[[t, Z]] of power series in the indeterminates (r, Z). 

1. This map is a morphism that associates to any f in K[Z] its Lie series defined 
by the formal power series J2 iefi T z 5 l f/i\. 

2. The derivation 5 is called the infinitesimal generator of e rS . 

3. The formal power series e rS Z are solutions of the vector field associated to S. 
These series form the formal flow of 5; this derivation induces an infinites- 
imal transformation from K x K m into K m that associates, under suitable 
condition of convergence, the evaluation (e rS Z)(V) to any parameter r inK 
and any initial point V in K m ; this map is the action of the flow e rS on K m . 

Example 7. Hence, the C(x)-morphism a T : x — > e T x could be defined by the 
exponential map ay := e rS where S denotes the derivation xd/dx acting on the 
field C(x). The set {a T r e C} is a one-parameter group of automorphisms. 

Lemma 8. Given two derivations d and 5, the Baker Campbell Hausdorff for- 
mula states that the relation e d e s — e s e d e^ d '^ holds. 

The next section presents some how some derivations could be expressed as 
translation in a suitable coordinates set. 

Some Algebraic Tools for Rectification of an Infinitesimal Generator. 

Principal element. Forthcoming manipulations are based on the existence of 
a special element that behaves as a time variable for considered derivation as 
shown by the following definition: 

Definition 9. An p element in a algebra A is principal for a derivation S acting 
on A if the relation Sp = 1 holds. 

In order to determine a principal element p of a derivation 8, one have to solve 
the following partial differential equation Sp — 1. As this is not a trivial task 
and as not every derivation has such an element, we are going in the sequel to 
restrict our manipulation to the following kind of principal elements: 

Lemma 10. Given a derivation S ofK[Z], if there exists an element g in K[Z], 

1. such that the relations Sg ^ and 5 2 g = hold, then the fraction p := g/Sg 

2. and a constant X of 8 such that the relation Sg = Xg holds, then for any 
constant c of this derivation, the transcendental element p := \og(cg)/X 

is a principal element p of 5. The element g is called the preprincipal element 
of S associated to p. 



Remark 11. As we adopt an algebraic standpoint in this note, this lemma re- 
quires to consider in the sequel the localization K[p, Z]^ of K[p, Z] at the mul- 
tiplicative closed set p := {(Spy \ i £ N} (resp. p :— {g 1 \ i 6 N}). In fact, given 
any canonical derivation d/dz of K[Z], there exists one, and only one, canonical 
derivation of K[p, Z], p s extending d/dz and such that the following usual re- 
lations dp/dz — dg/dz — (p/(8g) 2 )dSg/dz (resp. dp/dz = (1/ g)dg/dz) are well 
defined in K[p,Z}^ (for the sake of simplicity, we use the same notation for 
derivations acting on K[Z] and their extension to derivation acting on K[p, Z], p j). 
This shows that the Lie algebra DerKK[/9, Z]^ is well defined. 

Example 12. For any element h in K(Z), we consider the logarithm log/i i.e a 
transcendental field extension ¥L(Z, log h) and the associated derivation exten- 
sion such that Slogh = Sh/h. Hence, the derivation 5 := xd/dx acting on <C{x) 
has a unique extension to a derivation 8 acting on C(x,log(x)) such that the 
relation 8\og(x) = 1 holds. 

Construction of a Rectifying Coordinate Ring. Principal elements of a deriva- 
tion 5 allow to construct a rectifying field in which 5 acts as a simple translation. 

Lemma 13. Given a derivation S and one of its principal element p, let us 
define the following formal operator: 

^ p :=]T(-p)^. (12) 

As 5 is a derivation, this operator induces a homomorphism and the following 
exact sequence: 

- kcm s , p - K[p,Z] {p) ^4 K[[p,Z]] {p) - K{(] - (13) 

where the variables set ( denotes the set tts^Z of formal power series. 

Remark 14- To prove that the map ir$, p is a homomorphism, one can use the 
same argument then whose used in the proof stating the same property for the 
exponential map e s . By construction irg^p is equal to 0. Thus, the kernel ker^ p 
contains the ideal pK[p, Z]^ and is not trivial (see also Proposition 30). 

Using exact sequence (13), we could define a coordinate ring K[(] that is 
isomorphic to the quotient algebraK[p, Z]^/(ker 7T5 jP )K[p, Z]^ and a rectifying 
ring K[(, p] that is its finitely generated extension. To explain this terminology, 
first remark that the derivation 5 acting on K[p, Z}^ could be easily extended 
to a derivation acting on K[[p, Z]]^ and thus to K[p, £]. The following lemma 
states that the derivation S is rectified when we consider its action on K[p, (]: 

Lemma 15. With previously introduced notations, the following relations hold: 



Sp = 1, 5( = 0. 



(14) 



Sketch of proof . The first relation is the definition of a principal element. Ele- 
ments £ are defined by the series X^£n( — p) l 5 l Z/i\. By Leibniz' rule, we have: 

V*eN, ( 5((-p)^) = (-l)'(p^ + p- 1 ^), (15) 
and thus, derivation's linearity proves the last relations. 

The morphism ns, p induces a coordinates change allowing to express the deriva- 
tion S as a simple translation d/dp in this new coordinates set. We do not 
describe further this coordinates change, because we are just going to use some 
of its properties and not its exact formulation. Forthcoming considerations are 
based on the fact that the relations <5£ = imply that the relations e rS ( = ( 
hold. Thus, the morphism Tr$ tP maps the coordinate ring WL[Z] of the ambient 
space K m onto a coordinate ring invariant under the action of the flow e rS . 



2.3 Expanded Lie Point Symmetries and their Determining System 

Let us define now the derivations used in the sequel. 

Definition 16. Given a derivation 5, an algebraic system Sj and the associated 
derivation Dj, 8 is an infinitesimal generator of an expanded Lie point sym- 
metry of Sj if there exists a constant A of Dj such that the following relations 
hold: 

^ s Wte--dr-2^ f >-ter = xf « VieJ ' 

zez je.J J 

These relations form the determining system of Sj expanded Lie point symme- 
tries. 



(16) 
(17) 
(18) 
(19) 



Remark 17. Solution space structure — Derivations S satisfying (16) - (19) form 
a Lie sub-algebra of DerAK[Z] denoted by LieSym(Z'). Furthermore, if J\ is a 
subset of J 2 , the Lie algebra LieSym(Z',/ 1 ) is a sub-algebra of LieSym(Z , j 2 ). 

Remark 18. Considered Lie symmetries vs general Lie symmetries — The defi- 
nition 16 is designed for our algebraic purposes but is only a restriction of the 
general definition of Lie point symmetries (see [3]). In fact, remark that if the 
considered algebraic system Sj is 



a vector field (J = N), this definition reduces to the classical one of Lie point 
symmetries based on the Lie bracket i.e. [D,S] = XD; 

a purely algebraic system ( J = 0), this definition is more restrictive than the 
classical definition presented in Section 2.1 of [3]. In fact, let us consider the 
following system /i := x x 2 + y^ - 1, f 2 := x 2 2 + y 2 - 1, fz ■= x 2 y\ - y 2 x 1 
and the derivation S := x 2 d/dy\ — y 2 d/dx\ + x\d/dy 2 — y\d/dx 2 . The rela- 
tions 5 f\ = 2/ 3 , 5f 2 = —2/ 3 , (5/ 3 = fi — f 2 hold i.e. the derivation 8 leaves 
invariant the ideal spanned by fi,f 2 , fa and thus according to [3], 8 is the in- 
finitesimal generator of a one-parameter group of Lie symmetry (a family of 
morphism e rS parameterized by a constant r) that leaves the ideal spanned 
by {/1, f 2 , /s} invariant but not each of these polynomials (as shown by the 
relation e T& f\ = f\ cos 2 r + f 2 sin 2 r + 2/3 cos t sin r). This kind of deriva- 
tions are not taken into account in this note. This restriction is motivated 
by computational purposes i.e. the general definition of a Lie point symmetry 
for an algebraic system F = implies that we could — at least — determine 
in the quotient algebra K[Z]/FK[Z]; as this task is not in the complexity 
class considered in this note, we made a first restriction to the set of Lie sym- 
metry used in our work by only considering solution of system (16) - (19). 



Furthermore, there is little hope to solve the general partial differential prob- 
lem (16) -(19); thus, we restrict our solution space to affine infinitesimal gener- 
ators for which the associated determining equations form a linear system. 



3 Affine derivation and associated invariantization 



Definition 19. Let us denote by AffDerKK[Z] the following set of derivations: 

ls=^8(z)^- z \8( Zl ) :=b Zl +Y,a ZlZ2 z 2 , (b Zl , a ZlZ2 \ z 2 E Z) E K m+1 1 . 
I zez z 2 ez ) 

(20) 



Note 20. Notation — Given a derivation 8 in AffDcrKKfZ], we are going in the 
sequel to consider Z as a vector and use the following matricial notations: 



■As = (a ZlZ2 ) {ziiZ2)ez2 , B s = (b z ) zez , SZ = A S Z + Bs. (21) 



3.1 Determining System defining Affine Infinitesimal generators 



Lemma 21. For an affine infinitesimal generator 5 in KSDer^K[Z], the asso- 
ciated determining system (16) -(19) reduces to the following linear system: 



( 



df„ 





\ 

dh 

dz„. 



9fn 





(A 5 Z + B s ) - As 



V o ••• o / 



/ i \ 





V o J 



/l 

fn 

o 



(22) 



where H, e j is equal to 1 if the index i is in J and otherwise. 

Remark 22. Probabilistic resolution of determining system defining affine deriva- 
tion — The system (22) could be rewritten in a the more convenient matricial 
notation M(Z)K = where M(Z) is a m x (m + l)m matrix with coefficients 
in K[Z] and K is a vector whose (m + 1)™ coefficients are the coefficients of As 
and Bs- Affine derivations that are solution of this determining system, are given 
by the kernel of M(Z) in a field K. Kernel computation could be done by the fol- 
lowing probabilistic method. Indeterminates Z are specialized in matrix M(Z) 
to some random value in IK m in order to obtain a matrix Mi over the field IK; the 
resulting linear system M\K could be underdetermined and thus, several spe- 
cializations should be considered in order to obtain a linear system L t defined 
by M\K = ■ ■ ■ = M t K = 0. The rank r t of L l increases with i and the special- 
ization process could be stopped when r = r J+ i; the considered system L } could 
then be solved using a numerical method. The specialization set for which this 
process fails to find a correct solution is a zero-dimensional algebraic variety and 
thus, its probability of failure is low. 

However, there is an infinite way to choose a basis of the kernel computed 
above. But, one can use Lenstra, Lenstra and Lovasz' basis reduction algorithm 
in order to obtain a reduced basis in the sense that less variables are involved in 
each infinitesimal generators definition. 

To conclude, remark that some solutions of system (22) are spurious for 
our purposes since they describe the same flow and should be discarded. In 
fact, consider the problem £ defined by x = 9x: the base field of LieSym(Z') 
is the constant field of the derivation D := d/dt + 9xd/dx; thus, as d/dt is 
in LieSym(I7) and 9 is a constant of D, 9d/dt is another infinitesimal generators 
representing the same Lie symmetry then d/dt. These two derivations define the 
same Lie symmetry but are given by two different solutions of system (22). 



3.2 Principal Element Computation for Affine Derivation 

Lemma 23. Given a derivation 5 in AffDerKl&[Z] . If there exits a vector of 5's 
constants denoted by C :— (ci, . . . , c m ) such that the relations 



1. l CAg = and l CBg ^ hold, then the fraction t CZj t C{AsZ + Bg) 

2. l AgC = AC and 8\ = Q hold, then the element (log*C(Z + Bg/\))/\ 



is a principal element of 5. 

Sketch of proof . 1) Consider the polynomial CZ denoted by q. Using nota- 
tion (21), remark that Sg is equal to the linear combination t C(AgZ + Bg). 
Thus, the conditions on C given in the first item show that Sg is a constant 
different of and thus S 2 g is equal to 0. The first assertion of Lemma 10 is 
sufficient to conclude in that case. 2) Consider the polynomial t C(Z + Bg/X) de- 
noted by q. With the hypothesis of the second item, the element Sg is equal to Xg 
and thus, the second assertion of lemma 10 is satisfied and the transcendental 
element log(f?)/A is a principal element of S. 

Remark 24- Computational strategy — This lemma shows that in order to find 
a principal element for a derivation S in AffDerjKK[Z], one have first to check 
condition 1) and if it is not satisfied, one have to find an eigenvector of As. 



3.3 Flow of a Affine Derivation and Resulting Quotient Space 



Finding a coordinates change required to place a given derivation in rectified 
form is essentially the same problem as solving it in the first place. This could 
be easily done for affine derivation using the Jordan normal form as shown 
bellow. 

Remark 25. Jordan normal form — Given a m x m-matrix As associated to 
a derivation S, if its minimal polynomial p(£) is with p t — (£ — A t )" 1 

and Yl7=i a i = m i t nen there exists a change of coordinates V such that: 



As = V 



(■h 

J 2 

\0 ••• 



\ 





Jit: 



v- 



with J, := A, Id, 



/0 1 
00 1 

••• 
\0 ••• 



•0\ 



'•. 

1 
0/ 



(23) 



if Aj is different from and J t := A 4 Id Qj x Qj otherwise; the symbol Id Qi x Qj denotes 
the identity a t x ctj-matrix. This canonical form and Lemma 23 allow to compute 
preprincipal and associated principal element for any affine derivation 5. 

Hypotheses 1. From now, we suppose that the base field K is C in order to 
contain all eigenvalues of the matrix As and to define the quantities related to 
principal elements (exponentials and logarithms). 

Given an affine derivation S, one of its preprincipal element g and the associated 
principal element p, let us interpret geometrically the manipulation done in 
section 2.2. 



Flow associated to a derivation and induced equivalence classes. To do 

so and following Definition 6-3, we consider the application defined by the linear 
system of ordinary differential equations associated to our afhne derivation 8: 



& : 



(r,W) 



{e rS Z)(W) = exp(rA)W + ft exp((r - s)A s )B s ds. 



(24) 



this application could be computed numerically or using the following relations: 



01 T ' 



exp(rJ 4 ) = cxp(tA 2 ) 



K-i)! 







(25) 



1 / 



Note 26. Orbits of the flow — The image of K x W by constitutes an orbit of 
the flow e rS . This standpoint induces an equivalence relation ~ among the point 
of K m , with Vi being equivalent to V2 if these points lie in the same orbit of 
Let us denotes by Hs the set of equivalence classes a.k.a. the set of orbits. In the 
sequel, we suppose for the sake of conciseness that all the orbits have the same 
dimension i.e. we implicitly exclude from our statements the lower-dimensional 
orbits associated to the variety in K m defined by the ideal {S(z) = | z 6 Z}. 

All forthcoming manipulations rely on the following remark: 

Remark 27. Invariantization - Any object (algebraic relations, derivations, 
etc.) that is invariant under the action of the flow e rS will have a counterpart on 
the lower-dimensional variety H.$ whose representation — the invariantization of 
the considered object — completely characterize the original object. 

Note 28. As a first illustration, let us remark that any function f : K m — ► IK 
invariant for the flow e rS is invariant along its orbits and therefore there is a well- 
defined induced function f : Tl$ — > K; conversely, given a function f : TL$ — > K, 
there is an invariant function f : K m — + K defined by the relation f(V) := f (ft.) 
if V is in the orbit h. Hence, we obtain the following result: 

Lemma 29. There is a one-to-one correspondence between (polynomial) func- 
tions on K m invariant under the action of the flow e rS and arbitrary (polynomial) 
functions on Jig. 

In order to represent Tig, one can first find an algebraic representation of the 
orbits of the flow and then use a cross section of these orbits (see [5] and refer- 
ences therein for more details). This general approach is based on Grobner basis 
computation and treats general problems that exceed the scope of this note. 
Instead, we are going to use Lemmas 13 and 29 in order to give an algebraic 
description of Hs whose computation is based mainly on Jordan decomposition. 



Algebraic Representation of Hs- Consider the formal operator ns tP intro- 
duced in Lemma 13. Lemma 15 implies that the image of irs lP is invariant under 
the action the flow e rS . Thus, by describing the kernel of irs,p, we obtain an alge- 
braic description of functions on K m invariant under the action of the flow e rS . 
Lemma 29 shows that this description induces an algebraic representation of Us- 
The following lemma recapitulates these points when S is an affine derivation: 

Proposition 30. Given an affine derivation S, one of its preprincipal element g 
and the associated principal element p, one can define a homomorphism irs, p on 
an algebra of constants K[(] of S using the formal operator (12) as follow: 

0^q p K[Z] (p) ^K[Z] {p) ^K[Q^Q m] 

p(z) — p(C) 

where £ is equal to tts, p Z and q p is equal to l CZ if the preprincipal element g 
is defined by the case 1) in Lemma 23 and to t C(Z + Bs/ X) — 1 otherwise. The 
set of equivalence classes Hs could be identify with the hyperplane V(q p ) of di- 
mension m — 1 defined in K m by the linear form q p . Furthermore, 

1. the map Tr$ p induces a projection that associates to any point V in K m 
s.t. p(V) # 0, the point (tt s , p Z)(V) := £, eN {{- p{Z)Y {V Z) / i\){V) ofV(q p ). 

2. the points composing an orbit of e rS are projected to a single point in V{q p ). 

3. the orbit of e rS passing through a point in V(q p ) is projected on this point. 

Sketch of proof. Consider a principal element p of 8 and its defining preprin- 
cipal element g such that Sg = pg (rcsp. Sg ^ and 8 2 g = 0). Then, the rela- 
tion ns^ p g = ge~ (resp. Trg, P g — 0) holds and thus ng iP g is equal to 1 (resp. 0) 
(to be more precise dirs, P g/dz is equal to for all z in Z and thus ir$, p g is in K). 
If we denotes q — 1 (rcsp. g) by q pi the ideal g p K[Z] ( - p ^ is include in ker TTs,p- Fur- 
thermore, as we suppose that the flow e rS acts regularly, its orbits have the same 
dimension 1 and thus, the associated invariant coordinates ring is of dimen- 
sion m — 1. Hence, the quotient algebra K[Z]^/q p K[Z]^ allows to describe 
algebraically K[(}. Let us now, see some properties of the associated algebraic 
variety V(q p ). The first above assertion is just the definition of a projection 
on V(q p ). 

2) Remark that as the flow e rS is an homomorphism and relation e rS S = 5e rS 
holds, we have p{e T5 Z) l (b l e T& Z) = p{e rS Z) l {e T& : 5< Z) = e T& (p(Zy6 l Z), for all in- 
teger i. In order to show that the relation p(e rS Z) = e rS p(Z) holds, remark that 
as e rS is a homomorphism, g(e rS Z) is equal to e rS g(Z). If the principal element p 
is equal to log(g)/X (resp. g/Sg) with 5g = Xg (resp. 5g ^ and 5 2 g — 0) then 
the transcendental element log(e T(5 g) / X (resp. e rS (g/Sg)) is equal to log(pe TA )/A 
(resp. (g + T8g)/S(g + r)) and thus, p(e rS Z) is equal to p + r which is also equal 
to e rS p. Above relations show that (ns, p e T5 Z)(V) is equal to (e TS Trs lP Z)(V). As 
the flow leaves the image of tts iP invariant (i.e. e TS irg, p Z = -ks :P Z), this quantity 
is equal to (ns, p Z)(V). Hence, two points V and (e TS Z)(V) in the same orbit 
of e rS are projected onto the same point of V(q p ). 

3) If W is in V(q p ) then q p (W) is equal to and thus, p(W) is also equal to 0. 
In that case, by construction (tts vP Z)(W) is equal to W. 



The orbits of e rS cross the hyperplane V(q p ) transversally. In the sequel, we 
denote this hyperplane by Hs- Let us show now how works our reduction process. 



3.4 Invariantization and Parameterization 

Reduction of Algebraic Systems. Consider the variety V(F) in K m defined 
by the ideal spanned by F in K[Z]^ (for the sake of simplicity, we suppose 
that F is prime). 

Parameterization of an algebraic variety invariant under the action of e rS . If 5 
is an affine derivation such that the relation 5F = \F holds, it is an infinitesimal 
generator of a Lie symmetry e rS that leaves the variety V{F) invariant as shown 
by the following relations holding for all W in K m , 



V := (e TS Z)(W), F(V) := (F(e r6 Z))(W) = (e rS F(Z))(W) = e rX F{W). (27) 



As shown by Proposition 30, the hyperplane Hs is a linear cross-section of the 
orbit of e rS i.e. a variety that intersects these orbits in a single point. Further- 
more, there is a variety W induced by V(F) in Hs such that the variety V(F) 
is the image of K x W by the action of the flow e T<5 . Let us described now W. 

Invariantization of purely algebraic systems. The variety W is defined by the 
intersection V(Fmodq p ) n Hs (if <7 P (VV) and (F mod q p )(W) are equal to zero, 
then the relations F(W) = hold). As q p is linear, a description of W is obtained 
by a simple substitution in the equations describing V(F) (compare with the 
replacement invariant studied in [5]) as shown by the following example: 

Example 31. Let us consider the following purely algebraic system: 



S: {y-bf +a 2 =Z 2 /4, (x - a) 2 + b 2 = l 2 /4, x 2 + y 2 = I 2 . (28) 



Using results of Section 3.2, we determine its expanded affine Lie symmetries: 



As Hs a preprincipal element of Si, the solutions Z of system (28) are represented 
by the parameterization Z = e rSl Zi where Zi are the solutions of an invariant 
system obtained by the intersection of (28) with the hyperplane 1 = 1: 



As Si and 82 form an abelian Lie algebra, this last derivation is an infinitesimal 
generators of a Lie symmetry of 8\ ; it could be used to reduce further the sys- 
tem (30) (the Lie algebra spanned by 81 and 82 is abelian and thus solvable). In 
fact, the linear form a — lb — x is a preprincipal element of 82 associated to the 



Vz e {x, y, a, b, I}, z = e T zi 
h= 1, 




(30) 



principal element log(a — lb — x)/I. Using symmetry 5 2) solutions of system (30) 
could be represented as follow: 

f (y 2 - b 2 ) 2 + a 2 2 = 1/4, 
V*i G {x 1>yi ,ai,bi}, Zl = e^z 2 , ^ (1 + 7b 2 ) 2 + b 2 2 -1/4, (31) 
x 2 = a 2 - Ib 2 - 1, [ (a 2 - 7b 2 - l) 2 + y 2 2 = 1. 

In this particular example, the positive dimensional system (28) is represented 

— by a zero-dimensional algebraic system (31) that furnishes initial values Z 2 

— to an explicit linear differential system whose solutions are Zi = e r,52 Z 2 (this 
system associated to the derivation <5 2 is simple enough to be explicitly solved 
in closed form but in more complicated cases it could also be considered as 
a black box representation solved by purely numerical methods); 

— these values 7*i constitute an initial condition set of the linear differential 
system — induced by the derivation 5i — such that resulting solutions Z pa- 
rameterize the variety defined by system (31). 

We show now that the same type of results exists for differential systems. 
Reduction of Differential systems. 

Hypotheses 2. Restriction on Symmetries Specific to Differential Case — Lie 
symmetries of a given vector field D acting only on its state variables (i 7^ 0) 
could be used for Lie based integration but not for the previous reduction process. 
Thus, we suppose that D have an expanded Lie symmetry that acts at least on 
one of its parameter (Dz — 0) and that there is an associated principal elements 
such that the associated linear form q p satisfies the relation Dq p = 0. 

Invariantization of an infinitesimal generator. Given any derivation D acting 
on K[Z), the sequence (26) induces a derivation D acting on K[£] such that the 
relation tts. p o D = D o tts_ p holds. The exponentiation of D (see Definitions 6) 
induces a flow e rD on the hyperplane Hs that is the invariantization of the 
flow e rD acting on K m . Under above hypotheses, the flow e rD is just the restric- 
tion of e rD on Hs; in fact, as Dq p — the relation e rD q p = q p holds and thus, 
the flow e rD maps any point of Hs to another point of this hyperplane. The set 
of orbits of e rD in K m is projected onto the set of orbits of e rD in Hs- Let us 
see now the condition on 5 and D that allows to parameterize the set of orbits 
of e rD in K m by the set of orbits of e rD in Hs and the map e T<5 . 

Parameterization of vector field D invariant under the action of the flow e T<5 . 
If 5 is the infinitesimal generator of a symmetry of derivation D, according to 
Definition 16, the relation [D, 5} = \D holds. The Baker Campbell Hausdorff 
formula (Lemma 8) shows that the relation e TlS e ( 1+TlX ) T2D e~ TlS = e T2D holds. 
This implies that any orbit of e rD in K m is the image of an orbit of e rD in Hs by 
the flow e rS . Let us explicit all the process described above through an example. 



Example 32. Consider a FitzHugh Nagumo model (see § 7 in [1]): 

a = i> = c = d = 0, x = (x - x 3 /3 -y + d)c, y/dt = (x + a - by)/c. (32) 

The derivation 5 := d/dy + bd/da + d/dd is an infinitesimal generator of the 
following one-parameter group y — > y + A, a — > a + bX, d — > d + X that is com- 
posed of symmetries of the system (32). As the relation Sd — 1 holds, d is a 
(pre)principal element of S and the solutions Z := {x, y, a, b, c, d} of system (32) 
are described by the parameterizations Z = e dS Z where Z are solutions of a dif- 
ferential system on the hypcrplane V(d); hence, Z are given by the equations: 

y = y + d, a = a — bd, x = [x — x 3 /3 — y)c, y = (x + a — by)/c. (33) 
4 Conclusion 

In this note, we consider the computation of affine expanded Lie symmetries 
of a given algebraic system and show how this system could be rewrite in an 
invariant coordinates set for these symmetries in order to reduce the number 
of involved parameters. As this process is based on the computation of Jordan 
normal form and numerical linear algebra, its complexity is quasi-polynomial in 
input's size and likely polynomial for the great majority of practical cases. 

Extension of the reduction process to more general types of derivations. The 
manipulation presented in previous sections for affine derivations could be used 
for non-affine symmetries that occurs in practice as shown below. 

Example 33. Let us consider the following algebraic system: 

1 - x + x 2 y = 0, b - x 2 y = 0. (34) 
One can check that the following infinitesimal generator 

S := x 2 d/dx + (1 - 2xy)d/dy + x 2 d/db (35) 
is associated to the following one-parameter group of automorphisms: 

^■■^^ 2/-(Wra)(*A-l) 2 , b^b + x^, (36) 

that is a one-parameter group of Lie point symmetries of the system (34) . Remark 
that p := b/(x(x — b)) is a principal element of (35) (Sp = 1) and thus, all that 
we have done previously could be repeated i.e. the invariant coordinate set: 

TT SyP x = x-b, ir s , p y = (yx 2 - b)/ (b - x) 2 , ir StP b = 0, (37) 

allows to represent the solutions set (x, y) of (34) as follow 

x = x + b, y = {yx 2 +b)/(x + b) 2 V(x, y) s.t. 1 - x + x 2 y = 0, x 2 y = 0. (38) 



The results presented here could likely be extended for more general types of 
derivations but we do not know if the associated computations are feasible. 



Acknowledgments. The author is grateful to G. Renault, E. Schost, and M. 
SafeyElDin for many pleasant and useful discussions related to this note. 

References 

1. Murray, J.D.: Mathematical Biology. Volume 17 of Interdisciplinary Applied Math- 
ematics. Springer (2002) 

2. Khanin, R.: Dimensional Analysis in Computer Algebra. In Mourrain, B., ed.: 
Proceedings of the 2001 International Symposium on Symbolic and Algebraic Com- 
putation, London, Ontario, Canada, ACM, ACM press (2001) 201-208 

3. Olver, P.J.: Applications of Lie groups to differential equations. 2 edn. Volume 107 
of Graduate Texts in Mathematics. Springer Verlag (1993) 

4. Gatermann, K.: Computer algebra methods for equivariant dynamical systems. 
Number 1728 in Lecture Notes in Mathematics. Springer New York (2000) 

5. Hubert, E., Kogan, I.: Rational invariants of an algebraic group action. Construction 
and rewriting. Journal of Symbolic Computation (2006) (to appear). 

6. Burde, G.I.: Expanded Lie group transformations and similarity reductions of dif- 
ferential equations. In Nikitin, A.G, Boyko, V.M., Popovych, R.O., eds.: Symmetry 
in nonlinear mathematical physics Part I. Volume 43 of Proceedings of Institute of 
Mathematics of NAS of Ukraine., Kiev, Ukraine (2002) 93-101 



